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Abstract A novel Penning trap tower consisting of 
five compensated cylindrical Penning traps is developed 
for the PENTATRAP mass spectrometer at the Max- 
Planck-Institut fiir Kernphysik in Heidelberg, Germany. 
An analytical expression for the electrostatic potential 
inside the trap tower is derived to calculate standard 
Penning trap properties like the compensation of an- 
harmonicities and an orthogonal geometry of the trap 
electrodes. Since the PENTATRAP project described in 
the preceding article aims for ultra high-precision mass- 
ratio measurements of highly charged ions up to ura- 
nium, systematic effects for highly charged ions inside 
the trap tower are considered for the design process as 
well. Finally, a limit due to remaining anharmonic shifts 
at large amplitudes is estimated for the resulting geome- 
try, which is important for phase-sensitive measurements 
of the reduced cyclotron frequency of the ions. 



1 Introduction 

The mass of atoms and nuclei are of great interest in 
many fields of modern physics due to its inherent con- 
nection with the atomic and nuclear binding energies. 
The list of applications of precise mass or mass-ratio val- 
ues ranges from tests of nuclear and atomic structure, 
determination of fundamental constants and main fea- 
tures of the Standard Model for elementary particles to 
measurements related to neutrino physics [T]. This gen- 
eral need for very accurate mass values has triggered an 
ongoing technical development in mass spectrometers. 
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Penning traps have proven to be the ideal devices for 
the determination of atomic masses, capable of reach- 
ing extremely high accuracy down to 5m/m ~ 10~^^. 
Worldwide nearly 20 Penning trap mass spectrometers 
are in operation or under construction P]. 
With a destructive time-of-flight cyclotron frequency de- 
tection method, the most precise mass values reached so 
far for stable atoms have a relative accuracy of 5m /m ~ 
6 • 10"^*^ [S!,^]. Even for extremely short-lived species 
(down to ti/2 = 8.8 ms) or very low production rates of 
less than one atom per second in the case of heavy and 
superheavy ions, relative accuracies down to 5 -10"* have 
been shown [51|5]. For light and stable atoms and with a 
non-destructive image current detection method, results 
with an even more striking accuracy below the 10^^^ 
barrier have been achieved, allowing, e.g., for the most 
stringent direct test of Einstein's energy-mass-relation 
E = mc^ [7] and the most precise mass measurement of 
4He [g. 

This outstanding accuracy in Penning trap mass spec- 
trometry can be reached, since the measurement of the 
mass m of an ion with charge q is done via a determina- 
tion of the cyclotron frequency 



1 1 q 

— OJc = n 

ZTT ZTT m 



(1) 



where B is the magnetic field strength of the trap's mag- 
net. In Penning traps, the free cyclotron frequency Vc can 
be determined by application of the Brown-Gabrielse in- 
variance theorem [9] 

i^!^i^l + i^! + >^l, (2) 

where v^, and are the eigenfrequencies of the 
charged particle in the presence of the trap's magnetic 
and electrostatic field. Eq. ([2| is immune to small mis- 
alignments between the electrostatic axis of the trap and 
the magnetic field and also to small ellipticities of the 
trap electrodes. An exhaustive review of the physics of 
single charged particles in a Penning trap can be found 
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in [To]. 

The accuracy of measurements of atomic masses or mass- 
ratios is limited by a variety of systematic and statisti- 
cal effects. Examples are magnetic or electric trapping 
field imperfections and their temporal stability, eigenfre- 
quency shifts arising from the detection system, image 
charge effects or relativistic shifts of the ion's mass them- 
selves (see e.g. pnil2j ). The most crucial effect setting 
limits to the precision of the measurements is expected to 
be the temporal stability of the magnetic field. There are 
several concepts in current state-of-the-art experiments 
to minimize this influence. A very obvious approach has 
led to the development of an advanced stabilization sys- 
tem for superconducting magnets [13|, where relative 
drifts of 5B/B « 2 • lO^^Vh were reahzed. The tem- 
poral stability is even more crucial, since high-precision 
measurements are always performed in a relative way, 
which means that only frequency ratios between two ion 
species are measured. Therefore, not only the measure- 
ment itself, but also the loading and preparation of the 
two species has to be performed very quickly. Taking this 
to the limit, a very elegant method has been developed 
for low charged ions, where both ions are prepared on 
a common magnetron radius and measured at the same 
time [13] ■ Thus, temporal drifts and noise act on both 
ions in the same way and cancel out in the frequency 
ratio determination. 

For highly charged ions it is not possible to measure the 
mass ratio of two ions with high accuracy simultaneously 
in the same trap, due to the strong Coulomb interaction 
and the resulting systematic shift. Thus, the ions have 
to be measured either alternately in one trap, or simul- 
taneously in two different traps. Our project PENTA- 
TRAP [15], which aims for measurements of mass ra- 
tios of highly charged ions with an accuracy of < 10~^^, 
will enable both options. For this purpose, five identi- 
cal cylindrical Penning traps will be used to provide the 
possibility of flexible measurement schemes and fast ex- 
changes of ions between traps. The latter option, where 
the cyclotron frequencies of the two ions are measured 
simultaneously in adjacent traps is certainly favorable, 
since it suppresses statistical variations in the cyclotron 
frequency ratio due to temporal fluctuations of the B- 
field. Additionally, one trap can be used to store an ar- 
bitrary ion and permanently monitor its cyclotron fre- 
quency to detect temporal fluctuations of the magnetic 
field or to serve as a reference for the stabilization of 
the trap potential. The possible measurement schemes 
of this trap setup and ion candidates are described in 
detail in the preceding article by J. Repp et al. [15]. Be- 
sides temporal fluctuations of the magnetic field, which 
are taken care of by the measurement scheme, numerous 
other systematic effects have to be taken into account 
and intensive studies of the complete measurement sys- 
tem have to be performed. 

The subject of this paper is the design of the five-Penning- 
trap tower, which implies the calculation of the ideal or- 



thogonal geometry of the compensated traps. Further- 
more, systematic effects for highly charged ions inside 
the traps like image charge shifts and the Coulomb re- 
pulsion between two ions in adjacent traps are considered 
for the design process. Finally, machining imperfections 
are included and limits for measurements at large ampli- 
tudes are estimated, which is of special importance for 
the measurement of the reduced cyclotron frequency. 



2 Cylindrical Penning traps 

The natural geometry of a Penning trap utilizes hyper- 
boloidly shaped electrodes to generate a quadrupolar po- 
tential in the trap center [TU] . In our case cylindrical trap 
electrodes with open endcaps fTO will be used. This trap 
geometry enables an easy injection of the particles or ex- 
change between different traps. Additionally, cylindrical 
trap electrodes can be machined with higher precision 
and the resulting trap potential can be calculated ana- 
lytically, which is crucial for the design process. There- 
fore, this geometry has become very popular in recent 
years and is widely used at various Penning trap experi- 
ments like high-precision mass measurements, see e.g. [TJ 
l2| and references therein, (7-factor measurements of the 
bound electron and the proton [T7] as well as the free 
electron and the positron [T5] , cooling of highly charged 
ions [inj or the capture of anti-protons f^D] and for the 
production of anti-hydrogen [21.22 . 
If very high frequency resolution has to be achieved, it 
is essential that the electrostatic potential is harmonic 
along the trap axis to ensure that the axial frequency is 
independent of the energy of the ion. Therefore, correc- 
tion electrodes are added to compensate, e.g., for ma- 
chining imperfections. This was first done for hyper- 
boloidal traps in [23] and also applied to cylindrical traps 
in }16j . However, an additional correction voltage gener- 
ally changes the curvature of the harmonic potential and 
therefore the ion's axial frequency, which is very unde- 
sirable for practical use. Thus, a special geometry can 
be found where the harmonic part of the potential is 
independent of the applied correction voltage and stays 
constant while tuning out anharmonicities |16) . 



2. 1 Potential of a cylindrical trap tower 

An analytical expression for the potential of cylindrical 
trap geometries can be derived in a standard way from 
the Laplace-equation with azimuthal symmetry 

For a potential invariant under z t-^ —z transformation, 
the general solution is given by [23] 



(j){p,z)= / dfcA(fc)/o(fc/o) cos(fcz). 



(4) 
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Here, k results from the separation of variables and Iq is 
the modified Bessel function of zero order and first kind. 
The coefficient A{k) is defined by the boundary condi- 
tions, which are given by the geometry of the trap elec- 
trodes and the voltages applied to them. Unfortunately, 
this integral cannot be calculated analytically. Assum- 
ing closed and grounded endcaps, the potential has to 
vanish at the ends of the trap and the integral can be 
represented by an infinite sum, given by 



z)= ^ A(A;„)/o(fc„p)cos(fc„z), 



(5) 



where k is replaced by the discrete kn = Here, I 
is the total length of the trap tower. As discussed for 
a single open-endcap trap in [3S] , the suitability of this 
approximation for open-endcap traps depends on the ra- 
tio of the inner radius and the length of the grounded 
endcaps. It gives results valid to better than 1 % in all 
electrostatic coefficients if the length of the endcaps is 
three times the inner radius. In our case, this ratio is 
even higher which therefore leads to much more precise 
predictions about the electrostatic coefficients, as will be 
shown later. 

The coefficients A{kn) can be obtained by using the or- 
thogonality relation of the cosine function and integrat- 
ing over z along the inner surface of the electrodes at 
p — a. The potential at p = a at the small gaps between 
the electrodes can be approximated in axial direction by 
a linear interpolation between the values of the two ad- 
jacent electrodes. This leads to a negligible error if the 
length d of the gaps is much shorter than the length 
scales of the electrodes [Ml. In doing so, the following 
expression results for a trap with grounded endcaps and 
m = (2k + 1) electrodes: 



(j){p, Z) = ^ A{kn)Io (knp) cos (knZ) 
n—l,odd 



(6) 



A[krj 
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COS (fc„Z2i) " COS (fcnZ2i-l) 

Iq {k^a) 



(7) 



As shown in Fig. [T] Ui are the voltages applied to the 
electrodes, Z2i and Z2i-i specify the ends of the indi- 
vidual electrodes and a is the inner radius. In our case, 
every trap consists of five electrodes, in particular of a 
central ring electrode RE with length Ij. = 2zi, two cor- 
rection electrodes CE with length Ic ^ z^ — Z2 and two 
grounded endcaps EC with length 1^ — z^ — z^. 
The potential in Eq. ([6| allows one to deduce analyt- 
ical expressions for the static electric field coefficients. 
For hyperboloidal traps, extensive studies of these co- 
efficients and the influence on the ion's motion can be 
found in pilll27j . In [TF these results are transferred to 
cylindrical traps. In all of these works the potential is 




Fig. 1 Illustration of the a five-electrode cylindrical trap in- 
side a trap tower. The inner radius is given by a. The lengths 
of the trap electrodes are marked with Ir — 2zi for the central 
ring electrode RE, Ic — — for the correction electrodes 
CE and = — z^ for the endcaps EC. The length of the 
gaps separating the electrodes is given by d. The potentials 
at the inner surface of the electrodes are given by Ui with a 
linear approximation at the gaps between the electrodes. 



expanded in spherical coordinates in terms of Legendre 
polynomials. Additionally, in order to make the expan- 
sion coefficients dimensionless, the polynomials are mul- 
tiplied by powers of the small ratio r/do, where r is the 
radial coordinate and c?o a characteristic trap parameter 

In our work we follow another ansatz as shown in |26| . 
where the trapping potential is Taylor expanded in cylin- 
drical coordinates [p, cp, z] at the center of the trap. 



Hp, ^) = EE J\([^y. dp^dzo- 

j=0 i=0 W / P 



p^z^-^ (8) 



(0,0) 



oo ] 



(9) 



j=0 i=Q 



The coefficients Cij_ 



(0,0) 



Uo i'-U-t)'. dp^dzo- 

malized to the ring voltage U{). Due to the symmetry of 
reflection in z, all terms with odd order have to vanish 
to fulflU the boundary conditions. Therefore, both i and 
j have to be even numbers. In the following, we will use 
only the pure axial coefficients 



co,j 



1 iaJ>(p,z) 



j! 



dz3 



(10) 



(0,0) 



to analyze the potential. Due to W'^4'{p, z) = 0, the 
knowledge of the potential on the axis also determines 
the radial behavior. These Taylor coefficients are only 
functions of the geometry and the applied voltages and 
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can, therefore, directly be used to analyze the harmonic- 
ity of the trap geometry. By comparison of the coeffi- 
cients of same order j to the Legendre expansion used 
in [in] , it can be shown that the two kind of coefficients 
simply differ by a factor of l/{2dQ). Therefore, all results 
in |16j concerning the ion's motion can be used and are 
shortly summarized below in terms of the coefficients de- 



fined in Eq. ([10|). 

Certainly, the lowest order coefficients in Eq. ^ are the 
most important, since the axial amplitude of the ion is 
much smaller than the trap dimensions. For an ideal har- 
monic trap, the axial oscillation frequency is given by 



2C2 



(11) 



The higher order coefficients determine the anharmonic- 
ity of the trapping potential. They result in an amplitude 
dependent shift in the axial oscillation frequency. For the 
leading anharmonic coefficients C4 and Cq the shift in the 
axial frequency is given by 



C4 



5c6 E,\ 



4 cl qUo J qUn ' 



(12) 



where Ez is the axial excitation energy of the ion. In a 
five-electrode cylindrical trap, the potential can be tuned 
to minimize this anharmonic shift by adjusting the cor- 
rection voltage Ui = Uc (see Fig. [T]) . Since the potential 
in Eq. ^ is linear in the applied voltages, the Taylor 
coefficients of this potential can be written as 



Td 



J' 



(13) 



where T = Uc/Uq is the "tuning ratio" between the ring 
voltage Uq and the correction voltage Uc- The ej coef- 
ficients determine the part which is independent of the 
correction voltage, while the dj coefficients defined by 
dj = dcj/dT contain the influence of /7c- If T is chosen 
to be T = T\^^^Q — —64/^4, the leading anharmonic co- 
efficient C4 cancels out. Due to the dependence of C2 on 
Uc, tuning of the anharmonic coefficients by adjustment 
of Uc is accompanied with an undesirable change of C2 
and, therefore, also of the axial frequency in Eq. (11). 



To avoid this, a geometry has to be found in which d2 
and, therefore, the dependence of the axial frequency on 
the correction voltage is minimized. This trap geometry 
is called orthogonal. 



2.2 Ion-ion interaction 

Since we want to store at least three ions in adjacent 
traps in the final measurement process of our high-precision 
Penning trap mass spectrometer, an important issue is 
the Coulomb repulsion and the possibly resulting shift 
of the free cyclotron frequency. 

To analyze this effect, the additional ion- ion potential 
0ii(r, r') induced by a charge at r' has to be calculated at 



the position r of the ion-of-interest. Thereby, for our case 
the distance between two trap centers {Az ^ 10 mm) is 
much larger than the typical axial oscillation amplitude 
of the ions (^ 10 fim) and all dynamical aspects can be 
neglected. 

The analysis can be performend by expanding the ion- 
ion potential in a Taylor series with cylindrical coordi- 
nates [p,ip,z], in a way similar to the expansion of the 
trap potential in the previous section: 



0-f^'°''(p, 2;,p',z') 

00 n - 

V V t 

^ ^ m\(n 

n=0 m=0 ^ 



d"(l)tt{p,z,p',z') 



my. 



(0,0) 



= E E Sm.n-mP"^"""- (14) 
n— m— 

In contrast to the trap potential, for this interaction 
there is no symmetry with respect to the z-axis at the 
position of the ion-of-interest. Nevertheless, by inserting 



Eq. ( 14 ) into the Laplace equation it can be shown that 



the following relation of the Taylor coefficients holds for 
potentials with azimuthal symmetry: 



itaylor , 

Oii =. 90,0 +.91,0 



P + goA ■ 
- .91,1 • pz 



with gi.o 



.92,0 • P 

Oip 

- 0, .91,1 = and ^2,0 



50,2 



1 



r.9o,2- 



(15) 



(16) 



Hence, in second order this potential has exactly the 
same shape as the quadrupolar trap potential [TB]. In 
other words, in second order the ion-ion potential acts 
as a small offset on the trapping voltage, which exactly 
cancels in the determination of the free cyclotron fre- 
quency Vc [lUj. 

The remaining effect is a repulsion of the ions induced by 
30,1 resulting in a shift of the axial equilibrium position. 
The new minimum can be calculated from the superpo- 
sition of the ion-ion potential and the harmonic part of 
the trap potential 4>{r,z) — UaC2(z'^ — /9^/2) defined in 
Eq. ([9|. The minimum is given by 

.90,1 

"2(C/oC2+go,2)' 



zo = 



(17) 



where C2 = U0C2 can be determined by Eq. (11) for 



every q/m ratio. The coefficients go.i and 170,2 can be 
calculated by the potential of a point-charge in a con- 
ducting cylinder with inner radius a and length L given 
by HSj 



.(r,r') = 



neoa 



Jn 



(18) 



X sinh ( 



2 -^mn <-^77i-|_i (-^mn) Sinh ^ ^ • Xmn) 

sinh ' ^ 
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Here, Xmn is the n-th root of the Bessel function Jm of 
order m and 2>/< are the larger or smaller value of z 
and z'. 

Certainly the largest shift is expected for a shallow trap- 
ping potential corresponding to a large q/m ratio. On the 
other hand, the Coulomb repulsion in the trap tower is 
strongly suppressed compared to free space due to the 
shielding of the surrounding trap electrodes. 
Since the homogeneous part of our magnet is limited to 
120 mm, the maximum length of one trap is given by 
24 mm. It turns out that at this distance and for a large 
inner radius of a = 5 mm, this shift in the origin of oscil- 
lation is neglibigle for all experimental situations. Even 
in the extreme case of a proton as ion-of-interest and a 
charge of g = 81-1- in the adjacent trap, the shift is in 
the order of zq ^ 10^^^ m. Therefore, ion-ion interaction 
can be safely ignored for our trap tower. 

2.3 Image charge effects 

As pointed out in [^[50] for hyperboloidal Penning traps, 
the shift of the normal-mode frequencies in a Penning 
trap due to the electric field of image charges induced in 
the trap electrodes is not negligible for highly charged 
ions. In the case of cylindrical trap geometries, the re- 
sulting electric field and, thus, the force on the ion can be 
calculated analytically by solving the Laplace-equation 
as shown in [3T]. The shift of the radial frequencies is 
given as a function of the trap radius a 



q 



(19) 



whereas the axial frequency is not shifted, since the im- 
age charge electric field is translationally invariant in a 
cylindrical trap. The influence on the determination of 
the free cyclotron frequency can again be calculated by 
using the invariance theorem 

(20) 



AcOr 



ATTeQinaQUj'^ 

The resulting shift is shown in Fig.[2]for different masses. 
Note that for typical experimental conditions with aj+ 3> 
oj^ ^ a;_ , Eq. ( 20 ) is independent of the charge q of the 



ion. From this point of view it is convenient to have an 
inner radius as large as possible. At a = 5 mm, the shift 
in the free cyclotron frequency is still at ^ 5 • 10^^" for 
heavy ions like ^"^Pb. 

At charge states q > 20+, the relative shift per mass 
unit u is below 2.5 • 10"^^ for masses between A — 100 u 
and A = 250 u. Thus, if reference ion and ion-of-interest 
have nearly the same mass, the image charge effect on 
the ratio becomes negligible. As mentioned before, main 
interests of our project are neutrino related mass-ratio 
measurements or measurements of binding energies of 
electrons in highly charged ions. In this type of measure- 
ments, the mass differences of interesting candidates are 
typically well below 1 u. 




4 5 6 

a (mm) 

Fig. 2 Relative shift of the free cyclotron frequency due to 
image charges induced in the trap electrodes as a function of 
the inner radius a of the trap electrodes. Parameters used for 
the calculation are B — 7 T and ujz = 2n ■ 600 kHz. For the 
calculation, the q/m ratio was chosen sufficiently high such 
that Eq. ( |20[ ) becomes independent of the charge state q of 
the ions. 

3 Design of the PENTATRAP tower 

The design of our trap tower consisting of five identi- 
cal five-electrode traps is based on the potential given in 
Eq. (|6| . The Taylor coefficients of this potential were cal- 



culated using Eq. ( 10 1. Since the required order of Bessel 
functions n strongly increases with the total length and, 
therefore, the number of electrodes, all calculations were 
done with only three adjacent traps. For three traps 
the result of the calculations in the central trap are 
not affected by adding another electrode on the outside. 
Therefore, more electrodes are not needed for the calcu- 
lations. Furthermore, this shows that the approximation 
of the potential in Eq. ^ , which vanishes at the outer 



ends of the trap tower (see discussion in Sec. 2.1 ), has no 
influence on the coefficients in the center and is therefore 
justified for our calculations. 

For the real trap tower, the situation of two neighboring 
traps is common to the three inner traps and the results 
of the calculations are valid for them. But, certainly for 
the outer traps the performance will be slightly worse. 
On the other hand this is no issue to be worried about, 
since these traps are not dedicated for high-precision 
measurements (see [TBj). 

Due to the linear dependence of the coefficients on all 
applied voltages, the performance of the central trap cer- 
tainly depends to a small degree on the ring and correc- 
tion voltages in the adjacent traps. Therefore, all calcu- 
lations for the design of the trap were done with identical 
voltages in all three traps. This is the most common sit- 
uation for the planned measurements since mainly mass 
doublets, which need almost identical trapping voltages, 
are foreseen to be measured. 

To take advantage of the complete length of the homoge- 
neous region of our superconducting magnet of 120 mm, 
the length of Z = 24 mm per single trap was chosen. 
From the discussion of the image charge shifts and the 
ion-ion interaction (Sec. 2.2), an inner radius of a = 5 



mm was fixed as a compromise between these two ef- 
fects. The length of the gaps between the electrodes was 
set to d = 150 /im due to machining reasons. 
In addition to finding an orthogonal geometry where d2 
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vanishes, care was taken to keep the anharmonic shift in 



rL,=o-T^lc,=o has to 



Eq. ( 12 ) as low as possible in the final geometry. To this 
end, the trap should allow tuning out C4 and the next 
anharmonic coefficient cg simultaneously. Therefore, the 
difference in tuning ratios AT 
vanish. 

In order to find a suitable geometry, using numerical 
methods the following three equations were solved si- 
multaneously to find the unknown parameters Ir, Ic and 
T: 



C4{a,d,lr,lc,l,T) = 0, 
C6{a,d,lr,lc,l,T) = 0, 

^2(0, d, Ir, I, Ic) = 0. 



(21) 



The numerical results are summarized in Tab. [TJ where 
the geometrical lengths are rounded to 1 fim for machin- 
ing. In our case, the machining precision of the electrodes 
including a gold coating is limited to ±5 /im. The elec- 
trostatic coefficients are calculated for the ideal geome- 
try at the tuning ratio T = Tj^^^p. The errors resulting 
from limited machining precision are simply calculated 
by error propagation to estimate the worst case. In this 
respect, the contribution of the error of the gaps between 
the electrodes is approximately 10 times larger than the 
contribution from the error in any other dimensions of 
the trap. Note that the higher order coefficients still can 
be suppressed experimentally by finding a new tuning 
ratio which leads to C4 = in the real trap geometry. 
In Fig. |3j a technical drawing of the resulting assembled 
trap tower is shown. 



Table 1 Geometric and resulting electrostatic properties of 
a single trap within the PENTATRAP tower. The ideal ge- 
ometry is found by minimizing d2, C4 and ce simultaneously. 
For machining, the geometric values are rounded to 1 fj,m 
with respect to the ideal geometry. For an overview, the elec- 
trostatic coefficients up to 8 -th order are calculated for the 
ideal geometric values at the tuning ratio of the ideal geom- 
etry T = T|^^^(j. Their errors are estimated for a machining 
precision of ±5 /im. Note that the values of the higher order 
coefficients can still be suppressed experimentally by finding 
a new tuning ratio with respect to the real trap geometry. 



a (mm) 


5 


a (mmj 


0.15 


Ir (mm) 


1.457 


Ic (mm) 


3.932 


le (mm) 


7.040 




0.881 


C2 (l/mm^) 


(-1.496 ±0.007) • 10~^ 


d2 (l/mm'^) 


0± 1.241 • 10"* 


C4 (l/mm'*) 


0± 4.199- 10"'' 


di (l/mm'') 


(8.406 ±0.001) ■ 10""* 


cg (l/mm'') 


0± 1.892- 10"^ 


de {l/mm^) 


(3.579 ±0.019) ■ 10"^ 


cg (l/mm*) 


(1.672 ±0.077) ■ 10"^ 


ds (1/mm*) 


(-1.196 ±0.010) ■ lO"'' 



3.1 Orthogonality 

As a measure of the orthogonality of the trap the ratio 
d2/c2 can be taken, which is below 0.01 in the worst 
case. Therefore, following from Eq. (Ill and Eq. (13) a 



relative terms using Eq. (11) and Eq. (13): 



change of the tuning ratio shifts the axial frequency by 
only 



Av^ 1 d2 Hz 
— ^ = - — Vz = ±2.490 — — 
AT 2c2 mUrnt 



(22) 



at i>z = 600 kHz, where the error is again determined 
by machining imperfections of ±5 fim. In this context, 
the expression mUnit simply denotes a change of 10"'^ 
in the dimensionless tuning ratio T, which is a typical 
step size while tuning out anharmonicities. 
Since typical narrow-band detection bandwidths are in 
the range of 10 to 100 Hz, the ion's axial frequency 
stays well in this range while tuning out anharmonici- 
ties in mUnit-steps of the tuning ratio. Therefore, the 
ion's signal is always visible in the tuning process, which 
is of practical importance. Furthermore, for uncorrelated 
changes in the voltages applied to the ring and compen- 
sation electrodes, the effect of d2 can be expressed in 




(23) 



Due to this, stability requirements of the correction volt- 
age Uc are at least 2 orders of magnitude less than for 
the ring voltage Uq- 



3.2 Compensation 

In the PENTATRAP mass spectrometer we plan to use 
only cryogenic non-destructive ion detection. The axial 
frequency Vz will be measured with the dip method f 32j 
at cryogenic temperatures of 4 K. Due to the small par- 
ticle amplitudes in the order of 10 /xm, anharmonic shifts 
are of no concern for this method. The same is true for 
the measurement of the magnetron frequency , which 
is planned to be performed with the double-dip method 
explained in 33 . 

Anharmonic shifts are of more concern if detection meth- 
ods with excited particle amplitudes are used. For our 
setup this is the case for the measurement of the radial 
frequency z^+. This concern is even more important, since 
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EC CE RE CE EC 




120.6 mm 

Fig. 3 Technical drawing of the complete trap tower consisting of five identical, cylindrical and compensated Penning traps. 
The electrodes are labeled RE for the ring electrode, CE for the correction electrodes and EC for the grounded endcaps. The 
length scales of the electrodes can be found in Tab. [T] The electrodes are separated by sapphire rings, shown in white. 



the typical hierarchy of the three trap eigenfrequencies 
is iy+ ^ Vz ^ Therefore, the accuracy achieved by 
mass measurements mainly depends on the accuracy of 
the measurement of h'-^-. 

This frequency will be measured indirectly by a fast 
phase-sensitive detection method such as the so-called 
Pulse-aNd-Phase method (PNP) [33] or a similar novel 
method ^S]. In both cases, the phase of v+ is transferred 
to the axial mode by side-band coupling, where it is read 
out by detection of the signal peak. The frequency of the 
modified cyclotron mode is then given by the derivative 
of its phase with respect to evolution time. 
In general, the gain in measurement speed of phase- 
sensitive measurement methods compared to Fourier- 
limited frequency methods is approximately given by 
the factor 27r/cr(0), where a{4>) is the resolution of the 
phase measurement [36 . For low amplitudes, the phase 
resolution is limited by thermal phase jitter resulting 
from the noise background of the detector. Therefore, 
as a premise for fast measurements for both methods, a 
sufficient signal-to-noise (SNR) ratio has to be ensured 
by exciting the ion to amplitudes far above the thermal 
noise floor. In this context, anharmonic shifts can result 
in an upper limit for the excitation of the amplitude. 
This limit will be roughly estimated for our trap in the 
following paragraph. 

In most experimental situations it is convenient to tune 
out the leading anharmonic coefficient C4 by choosing 
In this case, the remaining anharmonic 



T 



C4=0- 



shift is mainly given by the next anharmonic coefficient 
cg. For the ideal geometry, both coefficients vanish for 
the same tuning ratio and the anharmonic shift given 
by Eq. ( 12 ) therefore also vanishes. For the real trap, 



this will not be the case anymore since these two co- 
efficients vanish at different tuning ratios. However, it 
will still be possible to compensate C4 only. To get an 
estimate of the remaining anharmonic shift in the real 
trap, a specific geometry has to be chosen since all coef- 
ficients Cj{a, d, Ir, Ic, I, T) and therefore the tuning ratio 
^lc4=o ~ T\^^^Q{a,d,lr,lc,l) are functions of the real 



geometry. 

As an example, we chose the rounded values for the trap 
geometry given in Tab.[l] In this case, the 6-th order co- 
efficient can be calculated to C6(r|g^^o) ~ ±1.892 -10"^ • 
1/mm^ with worst case mechanical tolerances included, 
and the difference in the tuning ratio of C4 and cg is given 
by AT EE Tl^^^o - = 22.9 • 10"^. As a result, 

the remaining anharmonic shift at low energies given by 
Eq. ^ ranges from d'^v^/ {dT dE^) ~ 3.5 mHz/(mUnit 
meV) to 1.4 mHz/(mUnit meV) at = 600 kHz for 
masses of 100 u to 240 u, respectively. As mentioned 
before, this will not be a problem for the low energy de- 
tection of Vz and . However, for the phase detection of 
v+ via the excited axial amplitude zq (see [34|[35]). the 
following limitation arises: 

Due to finite reproducibility of the excitation pulse and 
different starting conditions in the cooled axial mode, 
the final excited axial amplitude zq differs between two 
consecutive measurements. A typical value reported to 
us is a ~ 10 % deviation [35]. This in turn causes a shift 
of the frequencies of two consecutive measurements due 
to the amplitude dependent anharmonic effects given by 
Eq. (12). These frequency shifts limit the resolution of 
a phase measurement to cr(0) — Alo^ • Tmeas |36j . where 
Tmeas is the measurement time for the axial motion. 
With our calculated value of cg, this resolution normal- 
ized to the measurement time T^eas is shown in Fig. [4] 
For the measurement time, typically ^ 3 times the cool- 

ing time constant = ^ [32] of the excited motion 
is used. This constant depends on the charge q and mass 
m of the ion. Rq is the parallel resistance of the detec- 
tion system described in [TS] for our case. D^ff is the 
so-called effective electrode distance of the pickup elec- 
trode [39j , which is a measure of the coupling strength of 
the detection electronics to the ion. In our case I?e// is 
approximately 31.4 mm for the endcap electrode, where 
the detection pickup is planned to be. 
For charge states above q — 20-1- and our detection 
system, the expected cooling time constant can be es- 
timated to be below 300 ms even for masses up to 240 u. 
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Fig. 4 Phase resolution (t((j!)) due to anharmonic shifts in the 
axial osciUation frequency, normalized to the measurement 
time Tmeas- The shifts are caused by a 10 % difference in the 
amplitude zq of two consecutive phase measurements (see 
text). 



This implies, if a phase resolution of 10° and therefore a 
gain in measurement speed of approximately one order 
of magnitude compared to frequency measurements is 
aimed for, we are limited to 300 /im axial amplitude at 
a measurement time of Tmeas ^ 1 s. On the other hand, 
this is of no concern since the signal strength scales with 
q [55] and therefore, at g > 20+ amplitudes of 300 //m 
should ensure more than sufficient SNR. 



3.3 Simultaneous measurement in adjacent traps 

The main advantage of our five-Penning-trap setup is 
the possibility to measure the cyclotron frequencies of 
two ions of different species in adjacent traps at the same 
time, and hence suppress the effect of magnetic field vari- 
ation on the ratio. But, due to different charge-to-mass 
ratios there are different trapping voltages needed for 
the individual ions while the electrostatics of the traps 
was designed for identical voltages (Sec. [3]). Certainly, 
deviations from identical conditions can affect the per- 
formance due to the close distance and large openings 



of the traps. According to Eq. ( 11 ), the difference in the 
trapping voltages is determined by the difference in the 
q/m ratio between the ions in two traps. 
In many cases, e.g. the measurement of Q-values of /?- 
transitions, which are interesting for the determination 
of the neutrino mass, the corresponding mass difference 
is very small. E.g. in the transition of ^*^Re— > ^^^Os, the 
Q-value is only about 2.5 keV. Thus, for the same charge 
state of the ions, almost identical trapping voltages are 
needed which differ only by ^ 0.7 fiY in this case. 
The most critical case in this respect are measurements 
concerning binding energies of electrons in highly charged 
ions in order to test predictions of bound state QED the- 
ory. In this case the difference in the q/m ratio of two 
charge states is much larger. For the measurement of the 
mass-ratio between Pb*^+ and Pb*^+ at q/m — 0.39, the 
difference is about A{q/m) w 1.2 %. The trapping volt- 
ages at = 600 kHz differ by AU w 0.15 V. At this 



voltage difference between the traps, the anharmonic co- 
efficients in both traps are about twice as large compared 
to the case of identical trapping voltages. Resulting from 
the discussion of anharmonic effects in the previous sec- 
tion, this will not limit our measurement resolution for 
high charge states of the ions. 

4 Summary 

In this article the design of a novel five-Penning-trap 
tower for the mass spectrometer PENTATRAP [TS] was 
presented. The potential of a five-electrode cylindrical 
Penning trap inside a trap tower was calculated from 
the Laplace-equation and the electrostatic properties an- 
alyzed in terms of a Taylor expansion of the potential. 
Special attention was payed to the orthogonality of the 
calculated trap geometry as well as the possibility to 
tune out the leading two anharmonic coefficients of the 
electrostatic potential simultaneously. The experimen- 
tal limitations due to anharmonic effects at large am- 
plitudes were estimated, which is important for phase- 
sensitive measurement methods of the reduced cyclotron 
frequency. Furthermore, the systematic effects of highly 
charged ions in terms of image charge shifts of the free 
cyclotron frequency as well as Coulomb interaction of 
two ions in adjacent traps were estimated for the result- 
ing trap geometry given in Tab. [T] 
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